How Phase Transitions induce classical behaviour 
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We continue the analysis of the onset of classical behaviour in a scalar field after a continuous phase 
transition, in which the system-field, the long wavelength order parameter of the model, interacts 
with an environment, of its own short-wavelength modes and other fields, neutral and charged, with 
which it is expected to interact. We compute the decoherence time for the system-field modes from 
the master equation and directly from the decoherence functional (with identical results). In simple 
circumstances the order parameter field is classical by the time the transition is complete. 
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I. INTRODUCTION 

The standard big bang cosmological model of the early 
universe assumes a period of rapid cooling, giving a 
strong likelihood of phase transitions, at the grand uni- 
fied and electroweak scales 0] in particular. 

In this talk we describe how phase transitions natu- 
rally take us from a quantum to classical description of 
the universe. Metaphysics aside, cosmologists rely on the 
fact that the relevant fields obey classical equations from 
early times, since it is not possible to solve the quan- 
tum theory directly. Fortunately, we have reason to be- 
lieve that (continuous) transitions will move us rapidly 
to classical behaviour. Classical behaviour arises in the 
following way: 

• Classical correlations: By this is meant that the 
Wigner function(al) W[tt, (j)] peaks on classical 
phase-space trajectories, with a probabilistic inter- 
pretation. 

• Diagonalisation: By this is meant that the density 
matrix p{t) should become (approximately) diag- 
onal, in this case in a field basis. Alternatively, 
we can demand diagonalisation of the decoherence 
functional. In either case a probabilistic description 
(no quantum interference) is obtained. 

• Stochastic behaviour: The decoherence functional, 
which provides the diffusion (noise) to diagonalise 
the density matrix also supplies the dissipation that 
enables the fields to obey probabilistic stochastic 
equations, which evolve into classical equations. 

From the papers of Guth and Pi ^] onwards, it has been 
appreciated that unstable modes lead to classical corre- 
lations through squeezing. On the other hand, we under- 
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stand diagonalisation to be an almost inevitable conse- 
quence of tracing over the 'environment' of the 'system' 
modes. 

Continuous transitions supply both ingredients, from 
which the classical equations follow. Firstly, the field or- 
dering after such a transition is due to the growth in am- 
plitude of unstable long-wavelength modes, which arise 
automatically from unstable maxima in the potential. 
Secondly, the stable short-wavelength modes of the field, 
together with all the other fields with which it interacts, 
form an environment whose coarse-graining enforces di- 
agonalisation and makes the long-wavelength modes de- 
cohere. 

What matters are the time scales. An ideal situation, 
which we shall show is possible, is that the theory be- 
comes classical in the sense above, before the transition 
is complete. However, to quantify this is difficult be- 
cause, with fields, we are dealing with infinite degree of 
freedom systems. One of us (F.L) has shown elsewhere 
^ how classical correlations arise in quantum mechan- 
ical systems that mimic the field theory that we shall 
consider here, and we refer the reader to that paper for 
the role that classical correlations play. Our concern in 
this talk is, rather, with diagonalisation, determined both 
through the master equation for the evolution of the den- 
sity matrix and the decoherence functional, whose role is 
to describe consistent histories. Stochastic equations are 
then a corollary to this same diagonalisation. 

This talk builds upon earlier published work by us and 
Diego Mazzitelli j^^SHi together with our contributions 
to the proceedings of the 2001 meeting in Peyresq 0, Q 
and we refer the reader to this earlier work for much 
of the basic technical details. We restrict ourselves to 
fiat space-time. The extension to non-trivial metrics is 
straightforward in principle. See the recent work of Lom- 
bardo i9|, which complements this. The developments 
since the last proceedings are our greater understanding 
of the use of trial configurations and slower quenches , 
the parallel use of the decoherence functional to char- 
acterise decoherence, and the extension of the theory to 
include electromagnetism llC|. 
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II. EVOLUTION OF THE DENSITY MATRIX 

The evolution of a quantum field as it falls out of equi- 
librium at a transition is determined in large part by its 
behaviour at early times, before interactions have time 
to take effect. To be concrete, consider a real scalar 
field 4>{x), described by a Z2-symmetry breaking action 

{^^' > 0) 
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with symmetry breaking scale rj^ = 6/i^/A. On heating, 
this shows a continuous transition, with critical temper- 
ature = 2rj^ . If, by virtue of the expansion of the uni- 
verse the system is very rapidly cooled (quenched) from 
T > Tc to T < Tc, the initial stages of the transition can 
be described by a free field theory with inverted mass 
— /i^ < 0. The state of the field is initially concentrated 
on the local maximum of the potential, and spreads out 
with time. This description is valid for short times, until 
the field wave functional explores the ground states of 
the potential. 

The (p-tield ordering after the transition is due to 
the growth in amplitude of its unstable long-wavelength 
modes, which we term (j)^{x). For an instantaneous 
quench these have wave-number k < fj, for all time. 
Although the situation is more complicated for slower 
quenches, until the transition is complete there are al- 
ways unstable modes. As a complement to these, we 
anticipate that the stable short-wavelength modes of the 
field (/)>(x), where 



(x) = 0<(x) + (t>>{x), 



(2) 



will form an environment whose coar se-g raining makes 
the long- wavelength modes decohere In practice, 

the boundary between stable and unstable is not crucially 
important, provided there is time enough for the power in 
the field fluctuations to be firmly in the long-wavelength 
modes. This requires weak coupling A ^ 1. Of course, 
all the other fields with which (p interacts will contribute 
to its decoherence, but for the moment we ignore such 
fields (before returning to them in the last section). 
After splitting, the action ^ can be written as 

S[<j>] = ^[0<] + 5[0>] + ^i„t[0<,0>], (3) 

where the interaction term is dominated 0, 0] by its bi- 
quadratic term 



5'int [<j)< 



d'^x 4>\{x)4>l{x). 



(4) 



The total density matrix (for the system and bath 
fields) is defined by 



and we assume that, initially, the thermal system and its 
environment are not correlated. 

On tracing out the short- wavelength modes, the re- 
duced density matrix 



Pi-[(l)<,<f><,t] ^ J '^(t>>p[<P<,(t'>,(t'<,<f>>,t], (6) 

whose diagonalisation determines the onset of classical 
behaviour, evolves as 

Pr[t] = J d(f>+,J dcf,-, Mt,ti] p,[til (7) 

where Ji.[i,fi] is the evolution operator 

J,[i,ti]= f ^*I?0-exp{i^cG[0^><]}. 

(8) 

0<] is the coarse-grained effective action, of the 
closed time-path form 

ScG[^+,r<] - S[cp+] - S[r<] + 5S[(p+A<]- (9) 

All the information about the effect of the environment 
is encoded in (55* [0^, (^^] through the infiuence functional 
(or Feynman- Vernon functional [T^') 



F[0+,(/.-]=exp{z<55[(/.+ ,0-]}. 



(10) 



5S has a well defined diagrammatic expansion, of the 
form 

+\{{sU<PtAt])~[{su^tAt])f} 

-j (5'int [0< , 0>] 5'int [</>< Ay]) 
+« (5int [0< , 0> ] ) (5int [0< , 0>] ) 

+ \{SL[r<. - [{SUr<. <!>>])]'}■ (11) 

The quantum averages of the functionals of the fields are 
with respect to the free field action of the environment, 
defined as 



li P>['/'>i,'/'>i,io] 



^>f / , ^<^> / '^^> 



+ 



(5) 



xexp-{5o[(/)+]-5o[0;]}S[0+, 

To lowest order A^ diagrams are one-loop in the short 
wavelength modes. 



A. The Master Equation 

Once the reduced density matrix has become approxi- 
mately diagonal quantum interference has effectively dis- 
appeared and the density matrix permits a conventional 



probability interpretation. To see how the diagonalisa- 
tion of Pr occurs, we construct the master equation, which 
casts its evolution in differential form. As a first approx- 
imation, we make a saddle-point approximation for Jj. in 
Eq.®, 



In (|12|l is the solution to the equation of motion 



(12) 



SReScG 



= 0, 



and 



(13) 
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with boundary conditions 0^(io) = 

It is very difficult to solve this equation analytically. 
We exploit the fact that, even if the universe is com- 
pletely homogeneous prior to the transition then, after 
the transition, causality requires [l^ that it be inhomo- 
geneous because of the finite speed at which the order 
parameter fields can order themselves. This is in contra- 
distinction to the usual adiabatic analysis in which (for 
the continuous transition that interest us here) the cor- 
relation length diverges at the transition. 

Since the field cannot be homogeneous in either of 
its groundstates (f> = rj ot (j) = —rj there is an effective 
'domain' structure in which the domain boundaries are 
'walls' across which ip flips from one groundstate to the 
other. Further, these domains have a characteristic size 
^, where = nko labels the dominant momentum in 
the power of the (/)-field fluctuations as the unstable long- 
wavelength modes grow exponentially. For simplicity, we 
adopt a 'minisuperspace' approximation, in which we as- 
sume regular domains, enabling (/><ci(a;, s) to be written 
as 

0<ci(f, s) = /(s, i)<J>(x)$(y)$(z), (14) 
where $(0) = $(^) = 0, and 

^{x + = -*(a;)- 

f{s,t) satisfles /(0,t) — (/)<i and f{t,t) — 0<f. We write 
it as 

/(s, t) ^ (?!)<iWi(s, t) + (t><{U2{s, t). (15) 

In we made the simplest choice for $(a;), 

^{x) = cos k{)X. 

Extensions to include more Fourier modes are straight- 
forward in principle, but our work in was sufficient to 
show that the results only depend weakly on the details 
of the domain function <I>(x) for few Fourier modes. In 
the light of the more qualitative comments made here, 
we refer the reader again to for details. On the other 
hand, the Ui{s,t) are solutions of the mode equation 
for wavenumber kg during the quench, with boundary 




FIG. 1: The field profile 1141 in two dimensions, with lattice 
size ^. Dark areas represent (j) ^ rj, light areas (f) ~ —rj. The 
boundaries are domain walls, with profile given below. 



conditions Mi(0,t) = 1, ui{t,t) = and U2{0,t) = 0, 

U2{t,t) = 1. 

In order to obtain the master equation we must com- 
pute the final time derivative of the propagator J,-. After 
that, all the dependence on the initial field configura- 
tions (j>f^ (coming from the classical solutions 0^) must 
be eliminated. Assuming that the unstable growth has 
implemented diagonalisation before back-reaction is im- 
portant, Jr can be determined, approximately, from the 
free propagators as 

Jo[t,ii]= r"'2?0^ r"'l?0<exp{i[5o(0+)-^o(0<)]} 

(16) 

where 5*0 is the free-field action. This satisfies the general 
identities 



(17) 

which allow us to remove the initial field configurations 
if)^ , and obtain the master equation. 

Even with these simplifications the full equation is very 
complicated, but it is sufficient to calculate the correc- 
tion to the usual unitary evolution coming from the noise 
(diffusion) kernels (to be defined later). The result reads 



where D is the diffusion coefficient and 



(18) 



for the jinal field configurations (henceforth we drop the 
suffix). The ellipsis denotes other terms coming from the 
time derivative that do not contribute to the diffusive 
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effects. V is understood as the minimal volume inside 
which there are no coherent superpositions of macroscop- 
ically distinguishable states for the field. 



B. The diagonalisation of 

The effect of the diffusion coefficient in driving the di- 
agonaliation can be seen by considering the following ap- 
proximate solution to the master equation: 



Pr 



; ; t] cxp 



-VA^ [ ds D{ko,s) 
Jo 



(19) 

where is the solution of the unitary part of the master 
equation (i.e. without environment). The system will 
decohere when the non-diagonal elements of the reduced 
density matrix are much smaller than the diagonal ones. 

The decoherence time sets the scale after which we 
have a classical system-field configuration, and depends 
strongly on the properties of the environment. It satisfies 



1 « VA' 



tn 



ds D{kQ, s), 



(20) 



and corresponds to the time after which we are able to 
distinguish between two different field amplitudes, inside 
a given volume V. 

To terms up to order and one loop in the h expan- 
sion (we continue to work in units in which h — ks = 1), 
the influence action due to the biquadratic interaction be- 
tween system and environment has real imaginary parts 



ReSS = 



d\ I d% Aix)Kix,y)j:iy), (21) 



and 



Im55 = -i Jd^x I d^yA{x)N{x,y)A{y), 



(22) 



where K{x,y) = \X'^e{t - t')l-mG>\{x - y) is the dis- 
sipation kernel and N{x,y) = ^X^KeG-^^(x,y) is the 
noise (diffusion) kernel. G^_^_{x,y) is the thermal short- 
wavelength closed time-path correlator. The UV singu- 
lar parts of the loop diagrams are implicitly removed 
by renormalisation, leaving the finite temperature parts 
which are 0{T^). We also have defined 



(23) 



for final state modes. 

Explicit calculation shows that D{ko,t) is built from 
the diffusion kernel N as 



D{ko,t)^ / ds u{s,t) F{ko,s,t) (24) 
Jo 



where 



m(s, t) 



U2{s,t) 



U2{t,t) 

ui(t, t) 



ui(s,t) 



F(fco, s, t) is constructed from the spatial Fourier trans- 
forms of the overlap of the diffusion kernel with the 
field profiles ^{x)^{y)^{z). For the single mode <I>(x) = 
sin fcox 



A2 
64' 



F{ko,s,t) = '-A-[ReG^l{0;t~ s) 



-^ReG^l{2ko;t~ s) + ^ReG^l{2V2ko]t - s) 



1 



ReGlli2V3kQ-t~ s)]. 



(25) 



In the integrand of 1)24(1 u{s,t) is rapidly varying, driven 
by the unstable modes, and F{ko,s,t) is slowly varying. 
For long-wavelengths fco ^ M we have, approximately, 

F{ko,s,t)^0{N{ko^O;t-s)), 

whereby 

I?(fco,t) w i^(fco,0,i) [ dsu{s,t). (26) 
Jo 

That is, the diffusion coefficient factorises into the 
environmental term F, relatively insensitive to both 
wavenumber and time, and the rapidly growing integral 
that measures the classical growth of the unstable system 
modes that are ordered in the transition. 

To be specific, we restrict ourselves to the simplest 
case of an instantaneous quench from a temperature T — 
0{Tc)>Tc, for which 

sinhK(t-s)] sinh(tJos) 

"1 " ■ u/ — ' Ms,t) = , (27) 

smh(a;ot) smh(worj, 

where tjg = /i^ — fcp w /i^. It follows that 

M(s,t) = cosh2[wo(i - s)], (28) 

from whose end-point behaviour at s = of the integral 
(I26|l we find the even simpler result 

D(fco,t) - ^i-^F{ko,0,t) u{0,t) - (XTc/Awnf exp[2^i], 

(29) 

assuming fito ^ 1. The 0{T^) behaviour of F derives 
from the thermal short-wavelength modes. 

For more general quenches growth is more complicated 
than simple exponential behaviour but a similar separa- 
tion into fast and slow components applies. 

We have omitted a large amount of complicated techni- 
cal detail (see Q), to give such a simple final result. This 
suggests that we could have reached the same conclusion 
more directly. 

We now indicate how we can obtain the same results 
by demanding consistent histories of the 4> field. 



III. THE DECOHERENCE FUNCTIONAL 

The notion of consistent histories provides a parallel 
approach to classicality. Quantum evolution can be con- 
sidered as a coherent superposition of fine-grained histo- 
ries. If one defines the c-number field (t>{x) as specifying 
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a fine-grained history, the quantum amplitude for that 
history is ~ e*'^'"^' (we continue to work in units in 
which h — 1). 

In the quantum open system approach that we have 
adopted here, we are concerned with coarse-grained his- 
tories 



V<j) e 



iS\, 



(30) 



where a[(/)] is the filter function that defines the coarse- 
graining. 

From this we define the decoherence function for two 
coarse-grained histories as 



(31) 

I?[a+,Q; ] does not factorise because the histories (fx^ 
are not independent; they must assume identical values 
on a spacelike surface in the far future. Decoherence 
means physically that the different coarse-graining his- 
tories making up the full quantum evolution acquire in- 
dividual reality, and may therefore be assigned definite 
probabilities in the classical sense. 

A necessary and sufficient condition for the validity 
of the sum rules of probability theory (i.e. no quantum 
interference terms) is [T^ 



ReV 



0, 



(32) 



when ^ cT (although in most cases the stronger 
condition 'D\a^ , a~\ w holds 15]). Such histories are 
consistent Il6|. 

For our particular application, we wish to consider as 
a single coarse-grained history all those fine-grained ones 
where the full field remains close to a prescribed clas- 
sical field configuration i^ci- The filter function takes the 
form 



aci 



(33) 



In the general case, a[(j)\ is a smooth function (we exclude 
the case a[4>\ = const, where there is no coarse-graining 
at all). Using 



J(j) = I d'^xJ(x)(j)(x), 



(34) 



we may write the decoherence functional between two 
classical histories as 



X a+[J+]a-*[J-], (35) 

where 



jwij+,j- 



V(l)+V(l>~ 



„j(S[0+]-S[0-]+J+</.+ -J-0-) 



(36) 

is the closed-path-time generating functional. 

In principle, we can examine general classical solutions 
for their consistency but, in practice, it is simplest to re- 
strict ourselves to solutions of the form (|14|l . In that case, 
we have made a de facto separation into long and short- 
wavelength modes whereby, in a saddle-point approxi- 
mation over J. In this way, we can see that the above 
expression is formally equivalent to the definition of the 
influence functional (see Ref. Tl'l for details). Thus, we 
may write 



V( 



i) exp{iScG[ 



}• 



As a result. 



\V( 



exp{-Im55[0+ ^^i]} 



(37) 



(38) 



For the instantaneous quench of (|27() , using the late time 
behaviour 6^, ^ e'^'^ft'n, Ini^S'f^i, oi'Iil takes the form 



Im 5S ' 



1^ 



t 

ds ds'e^''' e^''''F{ko,s,s'). (39) 

"'0 



From this viewpoint adjacent histories become consis- 
tent at the time for which 



1 



dt Im 6S. 



(40) 



At this level, after performing the stationary phase ap- 
proximation, it is equivalent to evaluate the decoherence 
time scale from the master equation (through diffusion 
terms) or directly from the decoherence functional (or 
the influence functional). 



IV. THE DECOHERENCE TIME 

We have used the same terminology for the time t]j 
since, on inspection, (|40|l is identical to H20(l in deflning 
the onset of classical behaviour. As we noted, in practice 
the use of the decoherence functional looks to be less re- 
strictive than the master equation, and we hope to show 
this elsewhere. 

For the moment what is of interest is whether , based 
on linearisation of the model, occurs before backreaction 
sets in, to invalidate this assumption. When all the de- 
tails are taken into account, whether from H15|l or H27|l . 
to satisfies 



1 = O 



/V exp(4^iD) 



or, equivalently 

exp(4/i<D) = O 



(41) 



(42) 
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For the rapid quenches considered here, Hnearisation 
manifestly breaks down by the time t* , for which {(p'^)f ~ 
r?^, given by 



(43) 



The exponential factor, as always, arises from the growth 
of the unstable long- wavelength modes. The factor T^^ 
comes from the coth(/3u;/2) factor that encodes the initial 
Boltzmann distribution at temperature T >Tc. 

Our conservative choice is that the volume factor V 
is 0{^~^) since (the Compton wavelength) is the 
smallest scale at which we need to look. With this choice 
it follows that 



exp2(t* - to) 



O(05), (44) 



where 6 — 



b++ (b^)/2fi, and = 10+ - 0<|/2Ai. 
Within the volume V we do not discriminate between 
field amplitudes which differ by 0{fi), and therefore take 
(5 = 0(1). For (j) we note that, if were to equal t*, then 
02 = 0{1/X) = OiT^/n^) > 1, and in general > 1. As 
a result, if there are no large numerical factors, we have 



to < t* 



(45) 



and the density matrix has become diagonal before the 
transition is complete. Detailed calculation shows Q 
that there are no large factors [l^ . 

We already see a significant difference between the be- 
haviour for the case of a biquadratic interaction with an 
environment given by 120(1 and the more familiar linear 
interaction, adopted because it can be solvable (e.g. ^3)- 
This latter would have replaced A/zi^ just by S, incapable 
of inducing decoherence before the transition is complete. 
Although linear environments can be justified in quan- 
tum mechanics, in quantum field theory a purely linear 
environment corresponds to an inappropriate digonalisa- 
tion of the action. 

We note that, once the interaction strength is suffi- 
ciently weak for classical behaviour to appear before the 
transition is complete, this persists, however weak the 
coupling becomes. It remains the case that, the weaker 
the coupling, the longer it takes for the environment to 
decohere the system but, at the same time, the longer it 
takes for the transition to be completed, and the order- 
ing (|45|l remains the same. This is equally true for more 
general quenches provided the system remains approxi- 
mately Gaussian until the transition is complete. 



A. Back-reaction 

In both calculations for the decoherence time we have 
been obliged to assume that free-field behaviour explains 
the exponential growth of the long-wavelength modes. 
In reality, we are thinking of 0<f as describing the 



symmetry-broken phase, with magnitude 77, the symme- 
try breaking scale (if we normalise |$(x)| to be unity at 
its maxima). It can be shown |l9l | that, for an instan- 
taneous quench at least, nonlinear behaviour that stops 
the exponential growth only becomes important just be- 
fore t* . To see this, we adopt the Hartree approxima- 
tion, in which the equations of motion are linearised self- 
consistently. With a little work we find that the theory 
only ceases to behave like a free Gaussian theory with 
upside-down potential at a time t b , where 



t* -tB =^0{fi-^). 



(46) 



It follows that > in our ordering of scales 11. 



V. LATE-TIME BEHAVIOUR 

When is valid, we see that pr becomes diagonal 
before non-linear terms could be relevant. Although we 
haven't discussed it here, classical behaviour has been 
achieved before quantum effects can destroy the posi- 
tivity of the Wigner function Wi-, which is enforced by 
the unstable modes. Really, our to sets the time after 
which we have a classical probability distribution (posi- 
tive definite) even for times t > tB- The existence of the 
environment is crucial in doing this. 

This result also justifies in part the use of phenomeno- 
logical stochastic equations to describe the dynamical 
evolution of the system field, as we will now discuss. As 
it is well known [ill l20l | , one can regard the imaginary 
part of SS as coming from a noise source S.{x), with a 
Gaussian functional probability distribution. 

(47) 

where is a normalization factor. This enables us to 
write the imaginary part of the influence action as a func- 
tional integral over the Gaussian field ^(a;) 



Vax)P[^] exp 



d'^x A{x)^{x) 



exp <; - - / d^x I d^y A(x) N{x, y) A{y) 



;48) 



In consequence, the coarse-grained effective action can 
be rewritten as 

5cG[0t,0<] = --An J V^P[^] exp|z5eff[0t0<,^]|, 

(49) 

where 



The functional variation equation 



Aix)ax) 



0, 



(50) 
(51) 
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"semiclassical-Langevin" equation for the system-field 

MM 



>S'in 



dReScG[<f>t,^<,^2] 



'><=4>< 



(52) 



The evolution equation for the reduced Wigner func- 
tional Wr now becomes the Fokker-Planck counterpart 
to fK^ . 

Each part of the environment that we include leads 
to a further 'dissipative' term on the left hand side of 
(|52(l with a countervailing noise term on the right hand 
side. Although the 0<(/)5, and 0^0> terms were ignorable 
in the bounding of Id, in the Langevin equations they 
give further terms, with quadratic (fy^^3 noise and linear 
(additive) noise respectively. 

For times later than ts, neither perturbation theory 
nor more general non-Gaussian methods are valid. It is 
difficult to imagine an ab initio derivation of the dissi- 
pative and noise terms from the full quantum field the- 
ory. In this sense, a reasonable alternative is to analyze 
phenomenological stochastic equations numerically and 
check the robustness of the predictions against different 
choices of the dissipative kernels and of the type of noise. 
Hitherto, pure additive noise has been the basis for empir- 
ical stochastic equations in relativistic field theory that 
confirm Kibble's causal analysis However, recent 

numerical simulations with a more realistic mix of ad- 
ditive and multiplicative noise has shown that domain 
formation is unchanged ,22J . 



VI. FURTHER ENVIRONMENTS: NEUTRAL 
FIELDS 

Finally, it has to be said that taking only the short 
wavelength modes of the field as a one-loop system envi- 
ronment is not a robust approximation. This is particu- 
larly so for the Langevin equation [23| . We should be 
summing over hard thermal loops in the (/)-propagators. 
To be in proper control of the diffusion we need an envi- 
ronment that interacts with the system, without the sys- 
tem having a strong impact on the environment. This re- 
quires us to introduce further deconfining environments. 
We are helped in that, in the early universe, the order 
parameter field will interact with any field x for which 
there is no selection rule. Again, it is the biquadratic 
interactions that are the most important. 

The most simple additional environment is one of a 
large number TV 1 of weakly coupled scalar fields Xa, 
for which the action is extended to 



where S'[(/)] is as before, and 

N 



(53) 



N 
a=l 



(54) 



where > 0. For simplicity we take weak couplings 
A ~ 5a and comparable masses ma — /i- The effect of a 
large number of weakly interacting environmental fields 
is twofold. Firstly, the Xa fields reduce the critical tem- 
perature Tc and, in order that ~ x+Yl, g ^ 
must take A -I- ^ ^ 1- Secondly, the single x-loop 
contribution to the diffusion coefficient is the dominant 
X-field effect if, for order of magnitude estimates, we take 
identical pa = g/VN, whereby 1 ^ 1/VN ^ ~ A. 
With this choice the effect of the (/)-field on the Xa ther- 
mal masses is, relatively, 0{1/^/N) and can be ignored. 
We stress that this is not a Hartree or large-N approxi- 
mation of the typ e that, to date, has been the main way 
to proceed I25I l26| for a closed system. 

Provided the change in temperature is not too slow the 
exponential instabilities of the (p-field grow so fast that 
the field has populated the degenerate vacua well before 
the temperature has dropped to zero. Since the tempera- 
ture Tc has no particular significance for the environment 
field, for these early times we can keep the temperature 
of the environment fixed at = 0{Tc) (our calculations 
are only at the level of orders of magnitude). As before, 
we split the field a.s (p — (j>^+(j)y. The x-fields give an ad- 
ditional one-loop contribution to D{kQ,t) with the same 
u{s) but a G++ constructed from (all the modes of) the 
X-field. The separation of the diffusion coefficient due to 
X into fast and slow factors proceeds as before to give 
a term that is identical to H29(l (or (|39|l ) but for its 
prefactor. 

Diffusion effects are additive at the one-loop level, and 
the final effect is to replace A^ in (|1T|| by A^ -I- 5^ > A^ , 
while leaving (|43|l unchanged. Although the relationship 
between Tc and A has been uncoupled by the presence of 
the Xa, the relationship (|44l) persists, with an enhanced 
right hand side, requiring that (|45|l is even better satis- 
fied. 



VII. CHARGED FIELDS 

Given that the effect of further environmental fields is 
to increase the diffusion coefficient and speed up the onset 
of classical behaviour, additional fields interacting with 
the (p field seem superfluous. However, the symmetries of 
the universe seem to be local (gauge symmetries), rather 
than global, and we should take gauge fields into account. 
We conclude with some observations from our work in 
progress |lO| with local symmetry breaking. 

Local symmetry breaking is not possible for our real (f> 
field but, as a first step [3], it is not difficult to extend 
our model to that of a complex 0-field. At the level of 
0(2) global interactions with external fields and with its 
own short-wavelength modes, everything goes through 
essentially as before. The main difference is in the choice 
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of single degree of freedom configurations. Writing 

(l^ix) = -^{(pi{x) + i(j>2{x)), 

we assume that the (f)a behave independently until back- 
reaction is important. The simplest single-mode approx- 
imation to the long-wavelength system field is 



</'i,<ci(^,s) = /i(s,t)$(a;)$(y)$(z), 



(55) 



say, and 



h..<ci{x, s) = /2(s, t)<i>{x + a)$(y + + c), (56) 

for some non-zero a,b,c. fa{s,t) satisfies fa{0,t) — 4>a,i 
and fa{t, t) — (f>a,f- We write them as 



/a(s, t) = (l)a,iUi{s, t) + (t>aSU2{s, t) , 



(57) 



as before. Whereas the classical mode ()14|l of the real 
scalar described a regular array of domain walls, separa- 
tion ^, defined by the zeroes (/)<ci(x) = 0, the complex 
4><c\{x) describes a regular array of line zeroes (the in- 
tersections of 4>i<_c\{x) = = 4>2<c\{x)), which will evolve 
into global vortices after the transition. Although our as- 
sumption of a regular lattice of vortices is an extreme sim- 
plification, the production of vortices with typical sepa- 
ration ^ is as we would expect [T^ . 

In fact, to date we have not even been as sophisticated 
as (jSSjl and ifSBI) . but have just taken periodicity in a 
single direction j^]- This is sufficient to see that the sys- 
tem decoheres before the transition is complete, with an 
almost identical relation (|44|l . Wc assume that the insen- 
sitivity of the prefactor F{ko,s,t) to the regular lattice 
in both (|24|) and (|39l) is equally true here. This will be 
examined elsewhere (10|. 

Local U{1) symmetry breaking is most easily accom- 
modated by taking the 0-field to interact with other 
charged fields x through the local U{1) action 



in which S[(j), A^] — 



(58) 



and 



(59) 



S[A,„x]^ J d^x{{D^xrD''x + m\*x} 



(60) 



We have taken a single x held. The theory H58|l shows a 
phase transition, and we assume couplings are such as to 
make this transition continuous. 

For simplicity, let us just take x to be the environ- 
ment to the system field 0, which we do not separate 
into short and long-wavelength modes. On integrating 
out the X"fisld environment, the reduced density matrix 
Pr [0+ ,A+,(l)-,A-,t] evolves as 

Pr[t] = J d4>t j d^r JdA+ J dAr j,[t^ti] p,[U]. (61) 



(We have dropped the indices on for clarity). Yet 
again, we make a saddle-point approximation, 

Mtf,U] « eMiScG[c^tvA+,c^-^.A-,]), (62) 

where the coarse-grained action Scg has the form 

ScG[cl^+,A+,r,A-] = S[cf>+,A+]-S[cf,-,A-] 

+SS[<j)+,A+,cb-,A-]. (63) 

As before, SS encodes all the interactions between the 
environment and the system. In H62(l (/)^ is the solution 
to the equation of motion 



SReScG SReS, 



CG 



5A^ 







(64) 



subject to 0+ = and = A , with boundary con- 
ditions 0^(to) = 4^ E^nd (jy^i{t) = (j)f, and similarly for 

We stress that we are not tracing over the electromag- 
netic degrees of freedom, but determining the indirect 
effect of the x environment on the field, mediated by 
electromagnetism. 

Again, for simplicity, we assume an instantaneous 
quench. The diffusion is again driven by the unstable 
(j) modes that, approximately as 

(□ - x) = (□ - x) = (65) 

for times s <t*. This unstable scalar (pd is the source for 
the classical electromagnetic field, A^[(s,x), satisfying 

d''F,^{s, x) + (e^l^eiP + e2G++(0))^^,ci(s, x) 



+e- 



J dt J d^ylmA^,{s-t,x-y)A''^i{s,y) 



(66) 



= JmIS'X), 

in the Lorentz gauge, where — —ie(f>*^dj^(f)ci- 

In m, G++ix-y) = Tr[T{x{x)xHy)Pxm is the hot 
X propagator at temperature T. The G+-f (0) term is the 
X-loop thermal mass contribution to the field. 

We interpret (|66|l as being the start of an expansion 
with solution 

A',ci(s,x)= J d^ydtD^,{s,t;^-y)r{t,y), (67) 

where D^i,{s,t;x — y) is the thermal A^-field propaga- 
tor in the x heatbath. We have ignored the oscillatory 
solution of Af^ to the homogeneous equation, since this 
will not induce the exponentially growing diffusion that 
we need for rapid decohcrcnce. 

Just as for the other models considered earlier, when 
the transition is completed, there is a characteristic scale, 
the separation of the local vortices that express the frus- 
tration die to causal bounds. If we adopt a single char- 
acteristic scale before then, Im^S* now has two contribu- 
tions. We have already seen that the first, of the form 
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(|22|l . but from the x-loop, is sufficient to enforce deco- 
herence before the transition is complete, for acceptable 
parameters. We also have a contribution of the form 



ImSS 



-e 



d4y(AA)f^)ReA,,.(x,y)(AA)^^), 



(a)' 
(68) 



due to the electromagnetic field, where (AA)^ — 
A~j^ , derived from through H67|l . and 



y) 



f _d d_\ /_5 d_ 



X G++{x-w)G++{z-y) 



(69) 



This additional term to the diffusion function has 
derivative couplings. Having made a gauge choice, these 
give rise to explicit momenta factors fc^ in the generali- 
sation of F. Unlike the contributions to D that we have 
seen so far, which are largely insensitive to the momen- 
tum scale A:o, these contributions are strongly damped at 
large wavelength. In consequence, it is likely that they 
barely enhance the onset of classical behaviour but, given 
that the effect of the other environmental modes is to en- 
force classical behaviour so quickly, it hardly matters. We 
intend to give a fuller discussion of this elsewhere 0| . 



here before the transition is complete, under very general 
assumptions. An essential ingredient for rapid decoher- 
ence is nonlinear coupling to the environment, inevitable 
when that environment contains the short wavelength 
modes of the order parameter field. Had we only con- 
sidered linear coupling to the environment, as in |18| . 
for example (but an assumption that is ubiquitous in 
quantum mechanical models, from Brownian motion on- 
wards) decoherence would not have happened before the 
transition was complete, and we would not know how to 
proceed, although classical correlations would have oc- 
curred. For weak couplings further scalar environments 
with local interactions with the system field only make 
decoherence more rapid. However, it seems that, for the 
relevant case of a charged environment, also interacting 
indirectly through electromagnetic interactions, this in- 
direct contribution has little effect on a decoherence that 
is already effective. 
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